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^ Abstract 

^ . 

We discuss the construction of multi-caloron solutions with non-trivial holonomy, both 
as approximate superpositions and exact self-dual solutions. The charge k SU{n) 
> . moduli space can be described by kn constituent monopoles. Exact solutions help 
^ ! us to understand how these constituents can be seen as independent objects, which 
Q I seems not possible with the approximate superposition. An "impurity scattering" cal- 
I culation provides relatively simple expressions. Like at zero temperature an explicit 
■ parametrization requires solving a quadratic ADHM constraint, achieved here for a 
O ■ class of axially symmetric solutions. We will discuss the properties of these exact solu- 
tions in detail, but also demonstrate that interesting results can be obtained without 
explicitly solving for the constraint. 
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1 Introduction 

The last four years more understanding has been gained of the interplay between instan- 
tons and monopoles in non-Abelian gauge theories, based on the ability to construct 
exact caloron solutions, i.e. instantons at finite temperature for which Aq approaches 
a constant at spatial infinity [0, 0]. This last condition is best expressed by specifying 
the Polyakov loop to approach a constant value at infinity, also called the holonomy. 
It is the finite action that demands the field strength to go to zero at infinity, and 
guarantees the Polyakov loop to be independent of the direction in which we approach 
infinity. One can parametrize this Polyakov loop in terms of its eigenvalues exp(27ri/ij) 
and a gauge rotation g, such that 

Poo = hm Pexp( / Ao{x, t)dt) = exp(27rzdiag(/ii, /is, • • • , fJ'n))g, (1) 

which can be arranged such that Y17=if^i = 0, and /ii < /i2 < • • • < /Un < fJ'n+i, 
with fin+k = 1 + /Wfc- This is in the gauge where the gauge fields, assumed to be anti- 
hermitian matrices taking values in the algebra of SU{n), are periodic in the time 
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direction A^{x,t) = t + /?). In our conventions the field strength is given by 

We will find it convenient to construct the multi-caloron solutions from the ADHM- 
Nahm Fourier construction |1|, 0, ^, ^ based on taking instantons in R^, which peri- 
odically repeat in the time direction (see also Ref. p which uses directly the Nahm 
transformation ^ as the starting point for the charge 1 construction). To allow for 
non-trivial holonomy, the periodicity is only up to a constant gauge rotation (which 
is the holonomy). In this so-called algebraic gauge, all gauge field components van- 
ish at spatial infinity, and we may approximately superpose these calorons by simply 
adding the gauge fields. When each gauge field is periodic up to the same constant 
gauge transformation, the sum satisfies the same property. It should be noted that we 
are not allowed to add gauge fields with different holonomy; in an infinite volume the 
holonomy is fixed by the boundary condition. As to the topological charge, we recall 
that in the Atiyah-Drinfeld-Hitchin-Manin (ADHM) construction |^ it is supported 
by gauge singularities (the algebraic gauge is for this reason also called the singular 
gauge), as opposed to at infinity being a pure gauge with the gauge function having 
the appropriate winding number. To deal with the gauge singularity of one instanton, 
when adding the field of the others, one has to smoothly deform the gauge field of the 
latter to vanish near the gauge singularity. As long as the singularities are not too 
close, this can be done without a significant increase in the action. On the lattice this 
problem does not occur, when hiding singularities between the meshes of the lattice. 




Figure 1: Approximate superposition of two SU{2) charge 1 calorons with its pairs 
of equal mass constituents at x = (2,0,2), (2,0,8) and x = (8,0,2), (8,0,8). The 
logarithm of the action density is plotted as a function of x and z. The plot on 
the right shows one of the would-be Dirac strings, zooming in by a factor 40 on the 
transverse direction. 

Calorons, however, have an additional feature. When squeezed in the imaginary 
time direction (the size p becoming bigger than the period (3), they split in constituent 
monopoles with masses Sn'^Uj/P {uj = jij^i — jij). Outside the cores of these monopoles 
the gauge field becomes abelian. Ignoring the charged components, which decay expo- 
nentially outside the cores of the constituent monopoles, the field is described in terms 
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of self-dual Dirac monopoles. The singularity of a Dirac string can only be avoided by 
not neglecting the contributions coming from the charged components, even when far 
away from the constituents. For a single caloron we would not care, since the Dirac 
string is not seen in gauge invariant quantities. It does, however, involve a rather 
subtle interplay between the charged and neutral components of the gauge field in the 
vicinity of the would-be Dirac string [jl|. It is this subtle interplay that is disturbed 
when we add gauge fields of various calorons together. Unhke for the gauge singularity, 
the combined field will not diverge, but it shows a narrow and steep enhancement at 
the location of the would-be Dirac string as illustrated in Fig. |l|, where we added two 
SU (2) calorons. Also here one may shield the Dirac strings from these tails. But one 
always pays the price that the Dirac string no longer can be hidden, and carries energy. 
Let us stress again that these are genuine gauge invariant non-singular features in the 
configuration, even though they are of course a consequence of our particular way of 
constructing a superposition. 

A visible would-be Dirac string presents a formidable obstacle to considering the 
constituents as independent objects; they remember to which caloron they belong. 
Insisting the abelian field far from the constituents to be exactly additive under the 
approximate superposition leaves us little room for considering other possible superpo- 
sitions, apart from carefully fine-tuning the charged components of the configuration. 
Solving for the exact self-dual caloron solutions of higher charge we wish to show that a 
visible Dirac string is an artifact of the particular procedure to construct approximate 
caloron solutions. For the exact charge k caloron solutions we expect kn constituents 
and from the point of view of the parameter space these can be expected to be inde- 
pendent (as long the constituents do not get too close together). 




Figure 2: Comparing the logarithm of the action density (cutoff for log(S') below -3 
and above 7) as a function of x and z for the exact SU{2) solution (/i2 = \) with 
charge 2 (left) with the approximate superposition of two charge one calorons (middle) 
and the abelian solution based on Dirac monopoles (right), all on the same scale. The 
two pairs of constituent monopoles are located at x = (0,0,6.031), (0,0,2.031) and 
f= (0,0,-2.031), (0,0,-6.031). 

In this paper we develop the formalism and give exact solutions for a class of axially 
symmetric solutions. For SU{2) an example of such an exact solution with charge 2 
is presented in Fig. 0. We compare the exact solution (left), with the one obtained 
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by adding two charge 1 calorons, showing the would-be Dirac strings (middle) and 
with the exact abelian solution determined from (self-dual) Dirac monopoles placed 
at the location of the constituents (right). Indeed, the would-be Dirac strings are 
no longer visible for the exact non-abelian solution. This gives us good reasons to 
expect that moving away from the requirement of axial symmetry, the exact solutions 
will exhibit the constituent monopoles as independent objects. This is an important 
prerequisite for attempting to formulate the long distance features of QCD in terms of 
these monopoles. Non-linearities will remain important when constituents overlap, like 
for instantons at zero-temperature, as will be illustrated through suitable examples. 

Our results can also be used to extract information on multi-monopole solutions. In 
the charge 1 case it has been shown that sending one of the constituents to infinity 
one is left with exact monopole solutions. Perhaps somewhat surprisingly, it has been 
notoriously difficult to construct approximate superpositions of magnetic monopoles in 
such a way the interaction energy decreases with their separation. Within our approach 
this is due to the difficulty (particularly when sending some constituents to infinity) in 
keeping Dirac strings hidden. It would therefore still be important to find approximate 
superpositions that achieve this. Nevertheless, it shows how subtle it is to consider 
abelian fields as embedded in non-abelian ones. There is much more to 't Hooft's 
abelian projection 0] than meets the eye at first glance. 

The rest of the paper will be organized as follows. In Sect. 2 we formulate the 
Fourier analysis of the ADHM data for calorons of higher charge. This allows us in 
Sect. 2.1 to relate to the Nahm equations, providing the connection between the ADHM 
and Nahm data. The Nahm equation for charge k calorons expresses self-duality of 
dual SU{k) gauge fields on a circle, but with singularities. In Sect. 2.2 we introduce 
the "master" Green's function in terms of which the gauge field (Sect. 2.3) and the 
action density (Sect. 2.4) can be explicitly computed, relying heavily on some beautiful 
results 1^ derived in the context of the ADHM construction. 

This can all be achieved without explicitly solving the Nahm equation. In Sect. 3 
we will, however, solve it for some special cases with axial symmetry. There we also 
illustrate some interesting features that appear when constituents overlap, similar to 
what was observed at zero temperature [|TU|. 

Sect. 4 is devoted to the far-field limit (related to the high-temperature limit), 
where up to exponential corrections only the abelian component of the field survives. 
Again, expressions can be derived without explicitly solving the Nahm equation. We 
apply the formalism to the exact axially symmetric solutions in Sect. 4.3.1, where we 
prove that in the far-field region they are in a precise sense described by (self-dual) 
Dirac monopoles, which is conjectured to be the case for any exact solution. 

We end in Sect. 5 with a discussion of lattice results that have been obtained by 
various groups, some open questions and suggestions for future studies. 
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2 Construction 



In constructing the higher charge caloron solutions similar steps are followed as for 
charge 1. We distinguish two constructions, Sp{n) based on quaternions, and SU{n) 
based on complex spinors. With SU{2) = Sp{l), this implies that for SU{2) two 
constructions are possible. For k = 1 it was easily shown these are identical, but for 
k > 1 this is more complicated due to the quadratic ADHM constraint [Q] which can 
not be solved in all generality. Since the essential part of the construction that employs 
the impurity scattering technique is identical in both approaches, we will give a unified 
presentation. 

The SU (n) ADHM formalism for charge k instantons employs a k dimensional 
vector A = (Ai, . . . , A^), where A| is a two-component spinor in the n representation of 
SU (n) . Alternatively, A can be seen as a n x 2A; complex matrix. In addition one has 
four complex hermitian k x k matrices B^, combined into a 2k x 2k complex matrix 
B = using the unit quaternions = (I2, ir) and = (I2, —ir), where Tj are 

the Pauli matrices. With some abuse of notation, we often write B = a^j^B^. Together 
A and B constitute the {n + 2k) x 2k dimensional matrix A{x), to which is associated 
a complex (n + 2k) x n dimensional normalized zero mode vector v{x), 

A{x) = (^J^^^^ > B{x) = B -X, A\x)v{x) = 0, v\x)v{x) = 1„. (2) 

Here the quaternion x = x^cr^ denotes the position {a k x k unit matrix is implicit) 
and v{x) can be solved explicitly in terms of the ADHM data by 

u{x) = {B^ - x^)'^X\ (f){x) = ln + u\x)u{x), (3) 

As is an n X n positive hermitian matrix, its square root (f)^{x) is well-defined. 
The gauge field is given by 

Afj,{x) = v\x)d^v{x) = (f)~^ {x){u'^ {x)d^u{x))(j)^^ (x) + (f)^{x)dfj^(j)^^{x). (4) 

The 5*^(1) ADHM formalism for charge k instantons employs also a k dimensional 
vector A = (Ai, . . . , A^), where now Aj is a quaternion (4 real parameters). Again, A 
can be seen as an 2 x 2k complex matrix (but with 4k real, as opposed to complex, 
parameters). Now the four k x k matrices i?^ are required to be real and symmetric, 
still to be combined into a 2k x 2k complex matrix B = a^Bf^. The (2 + 2k) x 2k 
dimensional matrix A{x) is constructed as before. It is immediately obvious that now 
4>{x) is proportional to ctq, simplifying the expression for the gauge field to 

A^{x) = v^{x)d^v{x) = {u\x)d^u{x))/(j){x). (5) 
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For to be a self-dual connection, A(x) has to satisfy the quadratic ADHM 

constraint, which states that A^(a;)A(x) = B''{x)B{x) + A''"A (considered as a k x k 
complex quaternionic matrix) has to commute with the quaternions, or equivalently 

At(x)A(x) = ao®/-\ (6) 

defining as a hermitian (resp. symmetric) k x k Green's function. The self-duality 
follows by computing the curvature 

Ff,^ = 2(j)'^x)u'^{x)ri^^f^u{x)(j)'^x), (7) 

making essential use of the fact that /^^ commutes with the quaternions, and 1]^^^ = 
o'[fi(^u] being self-dual (f/^i, = (r^^a^j is anti-selfdual). The quadratic constraint can be 
formulated as Qm{A^{x)A{x)) = 0, where QmW = W — ^aQtT2W, and one obtains 

f]f,^ ® B^B^ + \Ta ® tr2(r„AU) = 0, (8) 

where tr2 is the spinorial trace. Note that this implies that tr2(TaA^) is traceless for 
a = 1,2,3. 

To count the number of instanton parameters we observe that the transformation 
A ATt, Bf, TB^T\ with T G U{k) (resp. T G 0{k)) leaves the gauge field and 
the ADHM constraint untouched. Taking this symmetry into account, one checks the 
dimension of the instanton moduli space to be Akn dimensional. We have Akn and 4fc^ 
real parameters from A and B^. The U{k) symmetry removes fc^ real parameters and 
finally the quadratic ADHM constraint gives real equations. On the other hand, 
for the quaternionic construction there are 4fc and A - \k{k + l) parameters from A and 
5^, of which \k{k — 1) are removed by the 0{k) symmetry, and the quadratic ADHM 
constraint gives 3 ■ \k{k — 1) equations. Global gauge transformations are realized by 
A — ^ g\, with g G SU{n) and are here included in the parameter count. For calorons 
with non-trivial holonomy the dimension of the gauge invariant parameter space is 
minimally reduced hy n — 1 (maximal symmetry breaking) and maximally by ra^ — 1 
(trivial holonomy). 

We note that for SU{2) the quadratic ADHM constraint and symmetry of the 
ADHM data for high charge differ considerably. For the caloron we have to deal in a 
sense with infinite topological charge, but finite within each (imaginary) time interval 
of length (3. This infinity is resolved by Fourier transformation, relating it to the Nahm 
formalism [Q, but the difference between the U{k) and 0{k) symmetries (the infinity 
is moved to making these gauge symmetries local) remains, as well as of course the 
nature of the Nahm data (the Fourier transformation of the ADHM data). Henceforth 
we put (3 = 1, which can always be achieved by a rescaling. 

Like for charge one the caloron with Polyakov loop Voo at infinity is built out 
of a periodic array of instantons, twisted by Voo- This is implemented in the ADHM 
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formalism by requiring (suppressing color and spinor indices, respectively quaternion 
indices) 

Upk+k+a{x + 1) = Upk+a{x)V;^^ (9) 

with p E Z. (the Fourier index), and a = 1, . . . , k (associated to the non-Abelian nature 
of the Nahm data). Using that (f)^2(^x + 1) = Poo^^^ (x)P^^, Eq. (H) leads to the 
required periodicity. Demanding 

^pk+k+a = 'Poo^pk+ay Bpk+a,qk+b = Bpk-k+a,qk-k+b + (^O^pq^ab-, (10) 

suitably implements Eq. (P) and is partially solved by imposing 

^pk+a '^ooCay Bpk-\-a,qk+b 

with A still to be determined to account for Eq. (H). It is useful to introduce the 

n projectors Pm on the mth eigenvalue of Voo, such that V^o = Z^m ^^'^^^"'-Pm. and 
\ = p2TTipti^ p /• 



2.1 Nahm setting 

We now perform the Fourier transformation to the Nahm setting [^], which casts B 
into a Weyl operator and A''' A into a singularity structure on S^, 



Y: Bpk+a,qk+b{x)e^^'^^-'^''^ = ^-i^-^b-\z'), 2m ^ e^-^^if = i'^^(z), 

p,q ■^'^'^ p 

E Aj,+.e2-(^^-^^')A,,+, = 5iz - z')Aabiz), ^ e-'""'"' Xpk+a = A.(z), (12) 



With B = a^B'^, we can write D{z) = a^D^{z) and A{z) = a^^Af^{z), where 

Df{z) = D^^z) - 2mx5''^ = S^^^o-^ - 2Tiix) + i"*(z), (13) 

Kb{z) = Yb{z- IJm)ClPmCb, Kiz) =Yj^(^~ l^m)PmCa- 

m m 

It should be noted that 5^ hermitian, implies that Ap_{z) is hermitian (as a. k x k 
matrix), whereas a real symmetric -B^ in addition implies A'jj^{z) = A^(— z). 

The 2x2 matrix A'^'' can always be decomposed as 

CPmCb = ^i<^os:!:-f-p:^). (m) 
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On the diagonal one can show, as for A; = 1, that = but in general the 

relation between and Sm is more complicated, 

qab qcd _i_ qad Qcb -^ab -*cd , -*ad -*cb ( ^ K\ 

m m ' m m rm rm ' rm rm ' V / 

Furthermore, tr2('fA^A) is traceless implies that J2a=i^m=i Pm — 0- Both these con- 
ditions are equally valid for the Sp{l) construction. But in the latter case, since 
through the Nahm equation pm determines the discontinuities in A{z), compatibility 
with Aj^^z) = A^{—z) requires pi + = 0, which will be verified below Eq. (|^). The 
Nahm equation is obtained by Fourier transforming the quadratic ADHM constraint 

D,{z)] T]^, = A7i'QmA{z), (16) 

or in more familiar form (as a k x k matrix equation) 

JzM^) + [M^).M^)] + ^,e,ke[Akiz),Mz)] = 2mY,6iz - p.m)pL (17) 
az ^ 

(A P2 PI pi-2 Pn-1 P{ Pi 

--_\ \ \ \ \ \ \--_ 

Z —pi p2 Ps Pn-2 Pn-1 Pn Pn+l=Pl + i- 

where the figure illustrates the jumps of Aj at each of the singularities (an overall factor 
of 27rz is not shown). The T symmetry in the ADHM construction translates into a U{k) 
gauge symmetry on 5*^, which allows one to set Aq = 27ii^Q to be constant, where is 
a hermitian matrix which can be made diagonal by a constant gauge transformation. 
Its trace part can be absorbed in xq. 



2.2 Green's function 

Central to the ADHM construction is the Green's function f^, which when Fourier 
transformed to f^''{z,z') = ^^^^ jp'=+«.9'=+''e2'^*(p^~'?^') is a solution to the differential 
equation 

+^T.Si^-Pm)Sm\w,z') = US{Z-Z'), (18) 

2m J 27T ^ J 

Si S2 S3 27rlfc Sn-2 Sn-1 Sn Si 

--J \ \---\r---\ — \ HH-z ^/ 

Z—pi P2 P3 Z Pn-2 Pn-1 Pn Pn+l=pl + ^ 

where the figure illustrates the jumps for the derivative of fx (itself being continuous) at 
each of the singularities, which now includes z = z' (an overall factor of 2n is omitted). 
With help of the impurity scattering formalism we will be able to express the solution 
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in a simple form, without assuming explicit knowledge of A{z). As a bonus this also 
gives a more transparent derivation for k = 1. The equation for the Green's function is 
valid for the SU{n) and Sp{l) formulations alike. In general the matrix is complex 
hermitian, but for 5*^(1) it is real symmetric, which implies that 

ff{z, z') = fl\z', z)\ and for Sp{l) only t\z, z') = /^(-/, -z). (19) 

For Sp{l) consistency requires, 5*1 — 5*2 = 0, which is to be compared with the condition 
we found in the previous section, pi + = 0. To demonstrate the validity of these 
relations we make use of the fact that (a = Ca^^^ with (j^ real (such that (I = (a)- 
Furthermore, we note that with Voo = exp(27rza; ■ f), /i2 = — Pi = uj = \uj\, whereas 
Pi = i(l2 — ci) ■ r) and P2 = |(l2 + ■ Introducing Q = lj ■ a = iou ■ r, which is an 
imaginary unit quaternion (i.e. Cl = —Q, QQ = co), we find with the help of Eq. (|1^ 

aoS:^-T.p^' = <%(ao-^(-l)™^])C,V, (20) 
= <aCb ([^{mM - «(-l)"'^[^fifT^]] + [^l,i(^u] - i(-l)™a{^fia^}]) 

for which we used that a^^^^cr,^] is a real quaternion, therefore given by icrotr2(a[^ncr,^]) = 
^o'Qtr2{firiu^) = (JoUJ-fjfj,iy We directly read off Sm and pm and verify that indeed 5*1 = 5*2 
and pi = —p2- We wrote the result so as to easily make contact with the k = 1 results 
(for which the 2nd term in Sm and the 1st term in pm do not contribute). 

In formulating the impurity scattering we note that the U (k) gauge transformation 

g{z) = exp (27ri(^o - Xolk)z) (21) 

turns D^i^) i'^to the ordinary derivative and conjugates all other objects in Eq. (p!8D, 
such that 

l^-^ + V{z;x)'^Mz,z')=A7rHk6{z-z'), (22) 
with fx{z, z') and V{z; x) given by 

fx{z,z')=g{z)Uz,z')g^{z'), V{z;x) ^ A7r^R\z;x) + 27rJ2Siz - f^m)Sm, (23) 

m 

where R{z; x) and Sm are defined by 

Rj{z; x) = Xjlk - ■^g{z)Aj{z)g\z), Sm = g{fJ'm)Sm9Kf^m)- (24) 

Periodicity is now only up to a gauge transformation. Note that Rj {z; x) is a hermitian 
k X k matrix, and that V{z; x) does not depend on xq. By combining and its 
derivative into a vector, 

Afe^')-(^{:t1)- (25) 
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we can turn the second order equation in to a first order equation 

Uz,z')=C5{z-z'), C= -47r2 



d 
dz 



Ik 
V{z; x) 



(26) 



Its solution is 

Mz, z') = W(z) \c{z') - 0{z' - z)W-'{z')C 





V{z; x) 







W{z) 



(27) 

where is a 2 x 2 matrix, whereas c is a two-component vector (hke C), all components 
being hermitian k x k matrices. The theta function takes care of the inhomogeneous 
part of the equation. However, we have to restrict to z' E [z — 1, z + 1] since the 
delta function is periodic. The solution can be extended beyond this range using the 
periodicity, which when imposed, as we will see, also determines c{z'). 

The solution for W{z) can be (formally) written as a path ordered exponential 
integral in the usual way. To show that this indeed makes sense, we should specify 
how to deal with the delta functions in V. From the definition of the path ordered 
exponential integral we find W{fim + 0) = TmW{fim — 0) with 



T 

J- rr 



Pexp 



fj.m+0 / g 



tk 




dz 



exp 





27iS„ 



( tk 





tk 



(28) 



which correctly reflects the matching conditions due to the "impurity" at z = fim, since 
fx{z, z') is continuous across z = /i^, whereas the derivative jumps with 2nSmfx{fJ'm, z') 
and both conditions can be summarized by fxif^m + 0, z') = T^fxil^m — 0, z'). Note 
that fx evolves with z as fx{z2,z') = W{z2)W~^{zi)fx{zi,z'), with 



W{z2)W-\zi) = Pexp [ r 





V{z;x) 



tk 




dz 



= W{Z2,ZI) 



(29) 



In particular, when zi = fim and Z2 = yUm+i this gives the "propagation" between two 
neighboring impurities and we can write 



Hm = W{flm+1 - 0, + 0) = Pexp 



Mm 





An^R^(z-x) 



tk 




dz 



(30) 



Neither nor if^ require us to specify a boundary condition for W{z). A change 
in boundary condition, however, affects c{z'). To avoid such ambiguities, we define 



W{zo)c{z') such that 

fx{z, z') = W{z, zo) \cM) - 0{z' - z)W-\z\ z^)C 



(31) 



We determine Czq{z') by scattering "around" the circle determined from the boundary 
conditions of fx- Using Eq. (|23|), the fact that fx{z,z') is strictly periodic and that 
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g{z+l) = gil)giz), one finds f^{z, z'+l) = f^{z, z')g^{l) and fxiz+1, z') = g{l)f^{z, z') 
(where the order of the matra multiplication, defined componentwise, is important). 
The latter condition can be used to fix Czo{z') over the range of one period, z' e [z, 1 + z] 
(the use of fx{z,z') requires z' e [z — l,z + 1], that of fx{z + l,z') further restricts 
the range for z' to z' > z). With the help of the periodicity properties of W it gives 

Czo{z') = (l2fc - g''{l)W{zo + 1, ;zo)) W'^^iz', zo)C. We thus find for arbitrary zq 

Uz, z') = W{z, ^o) {(l2fc - -^.o)"' - - z)^2k} W-\z\ z,)C, (32) 
where we introduced the "holonomy" J^^q 

J-,„^r(l)iy(zo + l,^o)=^7^(l)Pexp^^^ [y^^.^^ o'jrf^J- (33) 

The equation for is valid for z' G [2;, 2; + 1], but can be extended with the appropriate 
periodicity specified above. We note that ^^(1) plays the role of a cocycle, in the gauge 
where Aq{z) — 27iixolk is transformed to 0, with JF^^ the full circle "scattering" matrix. 
The 2;o-iiidependence of fx follows from 

= g^il)Wizo + l,z', + l)W{z'o + l,zo) = W{zo,z',)g\l)W{z'o + l,zo) 

= Wizo,z',)J^z>W-\zo,z',). (34) 

Putting things together we therefore find 

Uz,z') = g\z){lk,0)-Uz,z')g{z'), 
bl{z)Uz,z') = g\z)iO,lk)-Uz,z')g{z'), (35) 

satisfying all required conditions as can be checked explicitly. 

Interestingly, Eq. ( P^ implies J-'zg+i = 9{^)^zo9K^) should of course be the case. 
In the light of this we also note that (choosing zq = /i^ + 0) 

= TmHm~l ■ ■ ■ T2HiTig\l)HnTnHn~l ' ■ ' Tm+lHm, (36) 



" ^ ' t — r\ — ^ \ — 



using for the second identity that T^+n = g{^)Tmg\^) and Hm+n = .^(l)-f^m^^(l)- We 
will use these ingredients further on to relate to the earlier results for k = 1, where the 
positioning of ^^(1) is of course irrelevant. 
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2.3 Gauge field 

The central role of the Green's function f^{z,z') becomes clear when one appeals to 
the fact that it can be used to find the gauge field (working out Eq. (P), whereas its 
determinant gives a simple expression for the action density. This follows from the 
general ADHM construction P, and can be directly taken over for the caloron |jl], 
0, |TT|. For the gauge field one finds 

= \^'l\x)ff^^dA,{x)^^l\x) + \\r'l\x),d,<^^l\x)l (37) 
where 0(x) and <^j[x) are n x n matrices defined by 

0(x) = (l„ - A/,.At)"' , 0, = Aa,/,At. (38) 



To apply this to the caloron all we have to do is perform the Fourier transformation, 

m,m' m,m' 

(39) 

For the Sp{l) construction 0(x) is a real quaternion, and hence a multiple of do, 
after which the gauge field simplifies to 

A,{x) = \<P{x)ni,d,U^). (40) 

To simplify 0(x) and 0j(x) we use Eq. (|19|), together with /ii = — /i2, such that 
/f (/ii,/ii) = /^(/X2,/i2) and /f (/ii,/X2) = /^(/ii,/i2) = /^(/U2,/ii)*. We note that 
(j){x) and 4>j{x) involve the combinations Ca^^^Cb = CaCfiCb, which can be split in sym- 
metric and anti-symmetric combinations (cf. Eq. (PU])) 

Cacr.Cb = (TovUC^Cb+ClaajCb}. (41) 

No contributions from /"^(/ii,/i2) and /^''(Ai2, A^i) can appear in {/)(x) since these are 
symmetric in a and b, selecting from CaCoCfe fhe term proportional to do, but PiaoP2 
and P2<7oPi vanish. Therefore (cf. Eq. (^0|) ) 

0(X)-1 = ao- (/i2, /i2) El=l (PmCHaPm + (-l)"P„.r7^.Ca^C6 ^™) = (42) 



1 - /f (/i2, /i2) (eC + ■ V^uCaCi;]) = ^0 1 - 7r-iTr,(/.(/X2, /i2)5: 



Similarly we can simplify the expression for (f)j{x), which we split in a charged compo- 

bf{x) + (t)f 



nent and abelian, or neutral, component 0j(x) = (j)'p{x) + 0|'^°^(x) with 



=^ - /f (/ii, /i2)PlC{aa,a}P2 - h.C. (43) 



3 



12 



and 



/.abel / 



X 



/f (/^2, /i2) YL=1 {PmC{a^MPm + {-1)"^ PmVU^Q Pm) (44) 
= ■ T'/f (/^2,^2) [^tr2{iV ■ rCiaTjCb}) " i%uCaCb] = " ^TT-^Trfe (/^ (^2, /i2)P2) , 

where we used that itr2(u; ■ T({aTjCb}) = |tr2(C{6t^ • ^Caj^j) correctly identify pf", see 
Eq. (PUp. We may of course express (pix) and also in terms of the first impurity, 



TT iTrfc(/3.(/ii,pi)S'i 



and 



Aabel / 



XI 



lib ■ TTT ^Trfc(/^(/ii,pi)p{), as is 



</) \X} = (To 

easily verified. 

Like for = 1 we will show further on that /"''(/xi, /i2) decays exponentially, away 
from the cores of the constituent monopoles, where only the abelian component sur- 



vives, 



UJ ■ T 

2 



1 - TT ^Trfe(/^(/i2,/i2)5'2)l ?7f,^(9j. [vr ^Trfc(/a.(yU2,/U2)p2)l • (45) 



Note that for = 1 we are able to write A'f}"^^(x 



-2^ 



^ejr]l^^d^log(j){x), with 
e = P2/IP2I, niaking use of the fact that Sm = \Pm\, but that this is no longer true 
for higher charge, despite the fact that on the diagonal we still have S'^ = \Pm\- 
This seems to allow for the dipoles of k well-separated calorons to point in different 
directions. We will discuss this issue further when studying the limiting behavior far 
from the cores of the constituents, where the field becomes algebraic and constituent 
locations are readily identified. 



2.4 Action density 



Within the ADHM formalism the action density is given by 

- iTr„F2^(x) = -i9ja2log^(x), (46) 

where ip{x) equals 1/ det(/2^) after regularization to extract an irrelevant overall and 
for calorons divergent constant. We will be able to find a simple expression for iIj{x) 
at any k, generalizing the result for charge 1 calorons. We use that 9,^ log det(/^) = 
— Tr({c}y/~^}/a,) = -^Tr^D^fx), where in the last step we performed the Fourier trans- 
formation, and Tr includes an integration with respect to z. The case z/ = is treated 
separately, due to the discontinuity in the z derivative of fx{z, z') at z = z', which we 
regularize using point-splitting: 



lim 



dz Tri 



^fx{z + e, z') + ^fx{z - e, z') 
dz dz 



Airi / dz Tr 
Jo 



W~ 



[z, Zo) 



27rt J^'dz Tt {{l2k- J^zoY 





^l2fe 



W{z, Zo 

1 



(l2fc-^. 



J-2fc 



d 



Atc"^ dz 



TTk{fx{z,z)) 



(47) 
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The Tr without an index or hat indicates the full trace over the 2k x 2k matrix involved. 
To see how the total derivative term appears (not contributing to the integral due to 
the periodicity of Tvk[fx{,z, z)]) we use that 



^J^kiUz, z)) = -4vr2^Tr ( W-\z, z,) 





tk 



W{z,Zo) \{l2k-^zo)~^- St2k 



:-47r^Tr 



{w-\z, 





Ifc 



tk 
V{z; x) 



W{z,Zo) {l2k-J^zo) 



2k 



-At:'^Ty(w-^{z,Zo) 







li 



W{z, zo) [(l2fe-J^zo)" - st2k\ ,(48) 



with s = i (as for point-spitting, although one checks that the s dependent term 
actually vanishes). We note that J-'z^ depends on xo only through ^(1), and that 
dog{l) = 2TTig{l), such that doJ^zo — STri^F^g. With this we find 



do log det(/,) = -do log V', = det (ie-^'^° {l2k - ^.o)/^) , 



(49) 



which is independent of Zq. The factor i/\/2 in the argument of the determinant was 
inserted just so ip agrees with the definition introduced earlier for A; = 1. 

Next we compute dj log det (/c), 



1 

m Jo 



dz TrkiDiUz, z)) = -2 / dz IV^ (Rj(z; x)Uz, z)) = 



'' £ dz Tr (w-\z,zo) 
Trf(l2fc-^,o)"'-^.o r^'dz W-\z,Zo)d, 





Rj{z;x) 



W{z,zo) \{l2k-:Fj-'-sl 



^2k 



Ik 
V{z; x) 



Wiz,zo)), (50) 



where again s can take any value, but for convenience is best set to 1 here. Finally 
noting that only depends on x through R{z; x) and using that 



W-\z,Zo)d,W{z,Zo) - £^'°d^ W-\z,Zo)d, (^J. W{z,zo), 



(51) 



wc verify that 9^ log det (Z^:) = — S^log^^ for all /j, and k. It is amusing to note that 
this implies the remarkable formula 



iIV„F^2^(x) = -i9j9^1ogdet (l2k - ^^(l)Pexp 



Ik 
V{z;x) 



dz 



(52) 



even though explicit evaluation can be quite cumbersome. Not so for some special 
cases, including the single caloron k — 1, where R{z; x) is piecewise constant as we will 
discuss next. 
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3 Special cases 



— * — * 

Consider R{z; x) to be piecewise constant, for z G [nm, l^m+i] defined to be x — Ym (cf. 
Eq. (p4D), with Ym constant hermitian k x k matrices related to tlie Nahm potential 
by Aj{z) = 2nig'^{z)Y^g{z). In this case we can easily deal with the path ordered 
exponential integrals. To be specific, the "propagation" from n^. to yUm+i defined 
through Hm in Eq. (|30|), is given by 



^ ( COSh(27ri^m-Rm) i^irRm) ^ Smh{27[UrnRm) \ /ggx 

V27ri?„sinh(27ri/mi?m) cosh(27rz/mi?m) /' 

where z/^ = /i^+i - /"m (with ^m+n = 1 + Aim such that = 1) and Rm = 

{Rm ■ RmY^'^ (a hermitian k x k matrix). Since cosh(?/) and y"^^ smh(y) are both 
quadratic in y, actually no square root is involved in this expression for Hm- 



3.1 The known charge 1-caloron 

For charge 1 the Nahm equation (Eq. ([T7|)) has no commutator terms and R{z; x) is 
always piecewise constant. It reduces to the m*^ constituent radius for z G [/im,/im+i], 
Rm = rm = l^ml, whcrcas the prefactor at the impurity becomes Sm = I^Vi— "Tm-il = \pm\. 
With 



^ ^ / Tm \pm+i\ \ f cosh(27rz/mrm) sinh(27rz/mrm) 
~ ™ V r.m+1 ) \ smh{2-ni'mrm) cosh(27rz/„rm) 
a link to the earlier charge 1 results p| is established by noting that 



(54) 



'^'"-(27rw J)^'"+^^-(2vrr„ I) ' ^^^^ 

With the placing of ^^(1) irrelevant, and the possibility of absorbing in ^o, we 
therefore find 

= {2nrm o)'^^"""'^"^-i'^"^-2---AA---^m+i^m(27rr„ o) ' ^^^^ 

cf. Eq. ([361). In particular this shows that = — le"^'^*^" det(l2 — J^^^) agrees with 
the result found earlier, ip = itr(A„^„_i ■ ■ '-^i) — cos(27ra;o). The Green's functions 
can be shown to agree as well, using (I2 — ^^imY^ = (I2 — o'2^li^<^2) I det(l2 — ^/x„)- 

3.2 Exact axially symmetric solution 

Arranging R{z \ x) to be piecewise constant when k > 1 requires one to fulfill some con- 
straints. To solve the Nahm equation, Eq. (0), in terms of the Y^ = ^g{z)Aj{z)g'^ (z), 



15 



the commutator term should vanish. One way to achieve this, is by choosing Yjn = Ym(^. 
The Nahm equation relates the discontinuities of Aj{z) to pm, 

g\Pm){Yra-Ym-l)g{Pm)e = pm, (57) 

which imposes constraints on C,ai see Eq. (|14|). To seek a solution we choose all to 
be parallel in group space, = PaC {Pa a positive real number). This reduces the 
problem to k = 1, since (iPmCb = PaPbCPmC is proportional to C^P^C- For SU{2) 
this already solves the constraint, since for /c = 1 one has pi = — P2- For n > 2 it 
has been shown p that pm can take any value, provided I]m=i Pm = 0, in particular 
we may choose all pm to be proportional to e (by properly choosing Q. For k = 1 
it is convenient to parametrize ^PmC = {\Pm\ ~ ^' Pm)/(2vr) in terms of constituent 
locations, pm = ^Vm = ym — Vm-i- As in Sect. |2]^ we will take e = P2/|P2|- 

We obtain a larger class of for which the pm are parallel, by taking to be parallel 
up to a gauge rotation with an element of the unbroken subgroup f/(l)"^^ C SU{n) 
which leaves the holonomy unchanged, 

n 

(a = paexp{2niaa)(, = XI ^^^^rn, Tr„aa = 0. (58) 

m=l 

This leads to 

p"" = pVexp(27r^(«r - «a"^))Af™, ^^^f = p^l^ ■ AyJ\Aym\. (59) 

Note that for Sp{l) we verify that pi + = and S'l — 5*2 = (using al = — 
With = A?/me for all m (by definition Ay2 > 0) we may solve Eq. (p7|), 

" exp(27riy-a^, -(p, + sf)(eo^-eo)l) 



3 = 



' 2sin(7r[eo'-eo" 

(60) 

where the are arbitrary, m = 1, . . . , n and sj^ = i for j = 1, 2, . . . , m and = — § 
for j = m + 1, . . . ,n. The eigenvalues of these hermitian matrices determine the 
constituent locations, all lined-up along e. It should be noted that there is no reason 
to expect that all the Ym can be diagonalized simultaneously. We will come back to 
this in the following section. 

Returning to the simplest case of parallel gauge orientations, i.e. putting a{ = 0, 
we may take the limit ^ (related to vanishing time separations) to find 

n 

Y^ = i^a + VcpDSab + {Vm - yc)papb Vc = J2 ^rnVm- (61) 

m=l 

The {k = 1) "center of mass" coordinate yc can be freely chosen and the ^'^ play the role 
of "center of mass" of each constituent caloron. How exactly this is realized becomes 
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clear when we diagonalize Ym- Let us first consider Eq. ( |6TD for SU{2) and charge 2, 
k = n = 2, with = — ^2 = ^- Without loss of generality we choose i/c = 0, such that 
the two eigenvalues of Ym are given by = iy^p^ + (-l)^^^^ _^ i^2^p4 _^ y^^Ap'^, 

where = + p^ and Ap^ = Pi — P2- For large and positive ^ we find |/^^ = 
(— 1)-'^ + ymp'j + 0{^~^), representing two charge 1 calorons centered at ^ and — ^, with 
separations between their constituents monopoles given in terms of A1/2P1 2- We plot 
the constituent locations as a function of ^ in Fig. |^ for 7/2 = —yi = = h and pj = 2. 
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Figure 3: Constituent locations y^^ based on Eq. (^) (i.e. = and ^0 0) as 
a function of ^ = ^1 = —^2 for 1/2 = —yi = 1^1 = 1^2 = k Pi = P2 = 2. Dashed 
versus full lines distinguish the magnetic charge of the constituents. The dotted lines 
represent the four cases shown in Figs. H and |. 

Action density profiles are shown in Fig. ^ (left) for ^ = 3.5 and in Fig. ^ for 
^ = 1.6, 1.0, 0.5. From the dotted lines in Fig. | one reads off the associated constituent 
locations. Note that the magnetic moments of the two calorons are pointing in the same 
direction and that we can not freely interchange constituent monopole locations within 
our axially symmetric ansatz. However, when ^ is small it is more natural to interpret 
the configuration as a narrow caloron (i.e. instanton) with inverted magnetic moment 
in the background of a large caloron. This is the proper setting to understand the 
non-trivial time dependence for ^ = 0.5 illustrated in Fig. ^ (left). 



Figure 4: The action density (cutoff for log(S') below -3) as a function of x and z for 
the SU{2) solution with charge 2 (p2 = \, c^i = Q = Pj = 2) and increasing values 
of ^ = ^1 = -5 (see Fig g (left) for ^ = 3.5), with ^ = 1.6 (left), ^ = 1.0 (middle) and 
^ = 0.5 (right). Compare Fig. |^ for the corresponding constituent locations. 
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For ^ — ^ a singular caloron arises due to the fusion of two constituents (with 
opposite magnetic charge). This singularity is avoided when 7^ 0, which can be 
understood by observing that the eigenvalues of parametrize time-locations. If 
ai 7^ 0, with ,^0 and ^ made small, one will find two calorons (and their constituents) 
to be pushed far from each other. This can be understood as well, in terms of a short-to- 
long distance duality in the ADHM data for an instanton pair with non-parallel group 
orientation flOl, but can also be read off from the eigenvalues of defined in Eq. (|60|). 
As an example we take again SU{2) and charge 2, but now with ,^0 = ^0 = —^0 
a\ = — ttg = — «i = a2 = ct in general non-zero. For the case a = 1/8 (perpendicular 




Figure 5: The action density (cutoff for log(S') below -3) as a function of z and t 
(doubling the time-period) for the SU{2) solution with charge 2 (/i2 = |, ^ = |, a{ = 
0, pj = 2) and increasing values of ^o, with ^0 = (left), ^0 = 0.2 (middle) and ^0 = 0.25 
(right), for the top row all on the same scale, zooming in on the middle region on the 
bottom row (not to scale). See Fig ^ (right) for the case of ^0 = shown as a function 
of X and z. 



relative color orientations), pi = p2 = p and p2 



eigenvalues are y. 



VmP + 



(equal mass constituents), the 
+ i(A?/2)^p'^ sin~^(7r^o)) whereas for a = one 
In the light of this it is interesting 



finds ?/(^) = VmP' + i-iy^e + i(A|/2)V'cos-2(7reo) 
to observe, as shown in Fig. |^, that with a = and increasing ^0 the constituents are 
pushed out in the z direction as well. When ^0 ^ 0.5 the constituents would otherwise 
come close together through the periodicity in the time direction. Effectively these 
constituents thus have perpendicular color orientations (due to our choice of holonomy 
with P2 = j)- The transition from constituents separating in the time direction for ^0 
near to constituents separating in the z direction for ^0 near i occurs for p = 2 at 
approximately ^0 = 0.2. 

With a little imagination one detects the ring-shaped structure also observed in 
the case of instantons at zero temperature, see Fig. |^ (middle). A more direct analogy 
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of course occurs when two calorons (with p small, i.e. instantons with unresolved 
constituent monopoles) approach each other. We checked that for ^ = a = and 
^0 — *■ a singular caloron forms due to the overlap of two calorons with parallel gauge 
orientations, whereas for ^ \ the two calorons are pushed away (to infinity) in the z 
direction as is appropriate for the non-parallel group orientation due to the non-trivial 
holonomy. At an intermediate value (^o = 0.25 for pi = p2 = 0.1) one observes a small 
ring in the t-z plane. Choosing ^ large one may check that ±^o indeed gives the time 
location for each caloron. When, however, approaches \ they can no longer keep 
parallel gauge orientations due to the non-trivial holonomy. As noted before, this may 
be described by a solution with a 7^ and ^0 0. Computing the eigenvalues of Ym 
therefore allows one to easily predict the behavior of the exact solution. 

For charge 1 it had been shown that as soon as one of the constituents is far 
removed from the others the solution becomes static. For the "dimensional reduction" 
to take place at higher charge this is no longer sufficient. We have seen (generalization 
to SU{n) is straightforward) that any magnetically neutral cluster of constituents, 
when small with respect to /5, will behave like an instanton that is localized in time. 
For the special case with parallel group orientations, putting all = in Eq. (|6TD one 
would even be left with k — 1 singular instantons on top of one regular caloron, whose 
scale parameter is set by = Y!1=i pIi which can be understood from the fact that 
the matrix papb has rank 1. It does, however, give us one opportunity to go beyond the 
axial symmetry considered so far. When = we could solve the Nahm equation for 
parallel gauge orientations by = ymPaPb without insisting all the ym line-up. This 
still describes k — 1 singular instantons on top of one regular caloron, except that now 
the singular instantons can have arbitrary locations. 

Even though our ansatz to obtain exact solutions has been restrictive (as is clear 
from the axial symmetry), we stress that the solutions for the Nahm equation we found 
provide genuine multi-caloron solutions, which reveal kn isolated lumps for each of its 
constituent monopoles (with suitably chosen so the constituents do not overlap). 
This is not only illustrated in Fig. 0, but can also be understood analytically for any 
charge k and SU{n) as follows. We diagonalize each Ym with (in general different) 
similarity transformations Um, or Ym = Umdiag{yl^\ . . . ,yl{l^)Ul^. These bring Rm to 
the diagonal form Rm^^ = diag(r^'', . . . ,r^''), with r^-* = \x — ymS] the constituent 
radii, such that Hm = Ul^HmUm (Um acting componentwise) simplifies to (cf. Eq. (^)) 

~ ^/ cosh(27rz/^i?^-s) (27ri?^-s)-isinh(2vrz/„i?^-s)\ 
™ ~ V2vri?^^ssinh(27ri^^i?^^s) cosh(27ri/^i?^^g) )' ^ > 

The action density can now be explicitly expressed in terms of the constituent radii 

-iTr„F2^(x) = -19292 log^(x), ^ = det(^e-™o(l2fc - -^)/V2), 

T = exp(27r«(xolfc - ^o))UnHnfnHn-ifn-i " " " hMt,. (63) 
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cf. Eqs. ( p6| , ^ , ^ , |53D , where Tm = Ul^T.mUm~i {Um again acting componentwise). The 
size of the constituent monopoles is read off to be (2??//^)^^ (or (3{2iTh'm)~^ when /3 7^ 1), 
and one concludes that the action density will contain kn lumps for sufficiently well 
separated constituents. The figures were produced by computing the action density 
using precisely this method. 



4 Far-field limit 



The non-trivial value of the Polyakov loop at spatial infinity (holonomy) leads to a 
spontaneous breaking of the gauge symmetry, but without the need of introducing a 
Higgs field. One may view Aq as the Higgs field in the adjoint representation. This is 
one way to understand why constituent monopoles emerge. The best way to describe 
the caloron solutions, in case of well separated constituents, is by analyzing the field 
outside the cores of these constituents, where only the abelian field survives. Since in 
our case the asymptotic Polyakov loop value defines a global direction in color space the 
abelian generator in terms of which we can describe the so-called far-field configuration 
is fixed, giving rise to a global embedding in the full gauge group. Extrapolating the 
abelian fields back to inside the core of the constituents leads to Dirac monopoles. 
Such an extrapolation is well defined in terms of the high temperature limit, which 
makes the core of the constituents shrink to zero size and the field to become a smooth 
abelian gauge field everywhere except for the singularities of the Dirac monopoles. 
Thus we anticipate that in this limit the self-dual abelian field is still described by 
point like constituents, despite the fact that R{z,x) is no longer piecewise constant. 
Any "fuzziness" of the constituent location that may result from this, would be confined 
to the non-abelian core, and not visible from afar. 



4.1 Green's function 

Despite the somewhat formal expression for the Green's function fx{z,z'), one can 
extract information about the long-range fields from it. In the following we will show 
how to neglect the exponentially decaying fields in the cores of the monopole con- 
stituents, being left with the abelian components of fields which decay algebraically. 
We only need to consider the "bulk" contributions Hm, Eq. (0), which contain all 
the dependence on x. Our starting point is Eq. (^) restricted to the mth interval, 

|-^+4vr2^'(;^;x)|/™(^) = 0. (64) 

To distinguish between exponentially growing and decreasing contributions for this 
homogeneous equation we take as a basis for fm{z) functions f^{z) with the following 



20 



asymptotic behavior 

|f| oo : f^^{z) ^ exp {±27i\x\{z - firnjtk) y (65) 

relying on the fact that R'^{z; x) — > x^lfc. This prompts us to introduce on each interval 
the matrix valued functions R^{z) (we suppress the dependence on x) such that 



m 

fit 



±27r rR^{z)dz 



(66) 



from which it follows that -R^(-2) is a solution of the Riccati equation 

^^{^f ± ^™(^) = x). (67) 

We note that for |x| oo, R^{z) \x\ and that for piecewise constant R{z; x) both 
R'^{z) and R~Xz) are constant and equal to -R^, introduced in Eq. (|53D . 

We can write for z^z' G {fim, fJ'm+i) the "propagator" W{z,z') defined in Eq. (p9D 
in terms of /^(-z) as W{z,z') = Wm{z)W~^{z') with 

"^"^^'^ - [2nRl{z)mz) -2nR-{z)Uz) J ' ^'"^ 

Using that if^ = Wm{lJ'm+i)W~^ [fim) and f^{fim) = Ifc, we find by neglecting the ex- 
ponentially decreasing factors /~(/im+i) the required limiting behavior for Hm- Paying 
special attention to the ordering of the k x k matrices R^, observing that 

W;^\f.m) = (ArrRUf^^))-' QlR^tj X) ' ^"^ ^ '^(^trif^rn) + i?;;(^™)), 

(69) 

we find well outside the cores of the constituents 



The sparse nature of the matrices involved will be of considerable help to simplify the 
limiting behavior of the Green's function. A crucial ingredient is the combination 

f27rR-{fim) lk\f Ik 0\f Ik 0\ /27rS^ 

V 0j\2nSm ij l27ri?+_i(/x„) Oj \ 

for clarity summarizing the various ingredients in the following picture 

^=/^m-l R^_^{z),ji_^[z) f^m R^{z),f^{z) 
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This leads to the following far-field approximation for JF^^^ (cf. Eq. 



MiJ(^27r(i?++„_i(/i„+0 + ^^^+0 0j^'"l ^^^^ 
where Qm = Qm+n,m and 

Qm',m = /m'_l(/^m')(2-Rm'-l) "^^m'-l /m'-2 (A'm'-l) (2-Rm'-2) "^Sm'_2 

)(2i?„+i)-iS„+i/+(/i„+i)(47ri?„)-\ (73) 

One might have expected a factor 27rSm on the right, but this is contained in the 
remaining terms of Eq. ([7^). For example Trfc(^^^) = Tr^ (^5f^(l)^m(27rSm,)) • 

As we have seen in Eqs. ( P7| ) and the gauge field only requires us to know 
the Green's function at the impurities. Without loss of generality we may assume 
fim' > fJ-m and take zq = Hm + 0, such that (see Eqs. (^,^,^)) 

/x(Aim',/im) = -47r2(lfc,0) •Vr(/i„.,/i^ + 0)(l2fc-J^;,„)"^ f ) . (74) 



The matrix (l2fc — ^^j.^) has a 2 x 2 block structure, and one can verify that in general 



c d) ~ [ {b- ac-^d,)-^ (d-ca-^b)-^ 
Identifying the blocks, in the high temperature limit we find 
a = (l2fc - ^^lm)n Ifc - 27r^^(l)^„i?^(/i„;, 



(75) 



6= (l2fc- J-^„)l2 ^ -9\l)Qm, 

d={l2k-J'^,j22 ^ lfc-27r(i?+_i(/iJ + 5J^t(i)g^^ 



where we used that -R^+n- 1 (/^m+n )+Sm+n = g{X) {Rm- 1 (/im ) + 5"™ ) ^"^ ( 1 ) , cf . Eqs. (|4|,|6D . 
Evaluating the Green's function at the same impurities, /i^' = fJ^m, is simplified by the 
fact that iy(/im,/im) = l2k- This gives the following remarkably simple result in the 
far- field limit, 

/x(/in^,/i™) = -47r2(l2fe - J'^.Ju ^ {'^^^rn " 6^-'^(l))~'^ 27T (E^)"' . (77) 

For the Green's function evaluated at different impurities, fim 7^ fJ^m', we need to 
determine W{fim' — 0, /Xm + 0), for which we can follow the same method as for JF^^ 
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with Qrn',m as defined in Eq. ([73|). This leads to 



(1 



2k 



'^t^m ) 22 



-1 

Mm 7 12 



mi 



(79) 



which is exponentially suppressed since Gm',m grows as exp (27r|x|(/im' — yUm))- This 
cannot compensate for the decay of provided all fim are unequal, i.e. all con- 
stituents have a non-zero mass. Massless constituents have a so-called non-abelian 



cloud 12 , which has no abelian far-field limit. 



4.2 Total action 



To determine ^p in the expression for the action density, Eqs. (|46|J49|) , we need to 
compute det(l2fc — ^fi^)- Using Eq. ( [75D we find 

det {l2k - F^J = det f 5! Jj) = det f ^ = det(6) det{ab-^d - c) 



c d J \c — ab d d, 

^ det {g^{l)gm) det (g-'gil) - 27rS„) det (-27r^t(i)^^^s„) , (80) 

such that 

n 

^ ^ det(7r6;„S„) = 2-' n {det (/+(/i,n+i)) det(S„)/ det(2i?^)} . (81) 



m=l 



For |f| oo, f+inm+i) exp{2TrUm\x\lk) (see Eq. (|6|)) and iS„, i?^ ^ 
which implies that ip 2~'^nm=idet [exp(27rz/m|5^|lfc)] = 2^'^' exp(27rfc|a;|) (recall that 
J2m=i = 1)- Therefore, the action is given by 5* = — ^ / d'^x df^d'^ log ijj^x) = Sn'^k, 
as should be the case for a self-dual charge k solution. 



4.3 Gauge field 



Without the off-diagonal components of the Green's function contributing to the far- 
field region, the functions (j) and in Eqs (3^,39]) can be further simplified to 



(f){x 



-1 



mi l^'m^PmCa^jCb^'mi (^2) 



and only the abelian components of the gauge field survive. Particularly the case of 
Sp{l) discussed in Sect. p| is easy to deal with. Using Eqs. (|3],|0|,|2],|1^ we find 



[X 



[X 



UJ ■ (TTrfc[2S2 ^^(/i2)p2^^(/i2)]), 



u-f{l- Trfc[2S-^52])-^r7^,a.Tr,[2S-^^(/i2)p^2^T(/"2)]) 



(83) 
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where we recall (see Eq. ([71|)) that S2 = {Ri{fi2) + -R2^(/^2) + S2)- This is in perfect 
agreement with the earlier k = 1 results IQ. Note that the gauge rotation which 
relates fx{z,z') to fx{z,z') also relates Sm to Sm (see Eqs. (^,^)) and therefore does 
not appear in the final expression for (f){x). 

It is interesting to note that the dipole moment of the abelian gauge field is par- 
ticularly simple and does not require us to solve for R^{z), since lim|^|^oo ^2 = 2|x|lfc 
such that 

,lim A,{x)='-u.fv^^^d.^^^. (84) 

\x\^oo Z \X\ 

Hence the dipole moment p = iTrfc(p2) only involves (a, and we do except it allows for 
configurations with a vanishing dipole moment. For higher multipole moments, through 
R^{z), we need to deal with the full quadratic ADHM constraint, or equivalently 
with the Riccati and Nahm equations. Nevertheless, it is remarkable that in the high 
temperature limit the x dependence is restricted to Ri{fi2) and i?^(/i2) only. We 
would like to prove that each of its eigenvalues vanish at an isolated point, as one way 
to identify the 2k constituent locations. We will defer the study of this interesting 
issue, and its generalization to SU{n), to a future publication. 



4.3.1 Axially symmetric case 

The far-field approximation further simplifies when considering the axially symmetric 
solutions discussed in Sect. |3^ . We restrict ourselves here to Sp{l). Since R{z;x) is 
piecewise constant, the Riccati equation is trivial to solve. 



R^{z) = Rm= V(^lfc - eYm) ■ {xlk - eYj. (85) 
The square root involves a positive k x k matrix, and is well-defined. Due to the 

fact that Sm = g{l^m)Sm9Kl^m) = g{Hm)^ym " pm/ l^VmlgKl^m) (seC EqS. (^,^), the 

abelian component of the gauge field is of the simple form (e = P2/IP2I, see Sect. |273| ) 



AT i^) = --io ■ Te,%,dJog<P{x). (86) 

In the far field limit 0(x) — > (psix) (see Eq. (^)). Since S2 has rank 1, the matrix 
M = 2E2 ^5*2 has only one non-vanishing column with respect to a suitably chosen 
basis, which implies that 1 — Trfe(M) = det(ljfc — M). This allows us to write in the 
far-field region 

det(fii + Ri + g;) 
= det(fii + fi.-g.) ' 

from which we immediately read off the result [|T| for A; = 1, in which case it is easy to 
show that 0ff(x) = (r2 + e ■ r2) / {^i + e ■ ri), revealing A^{x) to be a linear superposition 
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of two oppositely charged self-dual Dirac monopoles. We would like 0ff(x) to similarly 
factorize for k > 1 in 2k Dirac monopoles, but since Ri, R2 and 5*2 in general do not 
commute, some care is required in demonstrating the factorization. We will rely on the 
fact that [Rm - e ■ Rm){Rm + e ■ Rm) = [Rm + e ■ Rm){Rm - e ■ Rm) = Rli- (e ■ Rm)^ = 
(x X e)^lfc and hence independent of m. Since by definition Ay2 > (see Sect, 
one finds that 

det ((i?i - e ■ ^1) + {R2 + e • ^2) 

0fT(a;) = 7 — 

det ((i?2 -e - R2) + {Ri + e-Ri 

This can now be reorganized according to 



6ff X = 



det ( 


\{Ri - e ■ Ri) + {R2 + e ■ R2)) 


1 {Ri + e-Ri) 


) det(i?2 + e-^2) 


det ((i?2 + e ■ i?2) ((^2 -e- R2) + {Ri + e- Ri)) 


) det(i?i + e-^i) 



_ det(i?2 + e-^2 
det(i?i + e- ^1)' 

after which we can separately diagonalize Ri and R2 to find 



B9) 



i Ti + e ■ Ti i Ti + r2 — \yi — y2 \ 

where y^e give the locations of the constituent monopoles, with y^ the eigenvalues of 
Ym and = x— y^^e (the index m distinguishing their charge). The second expression 
for the factorized version of 0ff(x) uses the fact that the constituent locations can be 
ordered according to 

,«<,«<,r<,f <...<, f)<,f. (91) 

It should be noted that this prevents the constituents to pass each other while varying 
the parameters for the axially symmetric configuration, see also Fig. ^ and the discus- 
sion in Sect. p.2[ We need to go beyond this simple axially symmetric configuration to 
allow for the constituents to rearrange themselves more freely. 



5 Discussion 



We have presented the general formalism for finding exact instanton solutions at finite 
temperature (calorons) with any non-trivial holonomy and topological charge. In an 
infinite volume holonomy and charge are fixed. The solution is described by 4kn 
parameters of which 3kn give the spatial locations of the kn constituent monopoles. 
The remaining parameters are given by k time locations and {n — l)k phases associated 
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to gauge rotations in the subgroup that leaves the holonomy unchanged. Of these, n — 1 
can be considered as global gauge rotations. The dimension of the moduli space, i.e. 
the number of gauge invariant parameters, is therefore equal to Akn — [n — 1). Our 
subset of axially symmetric solutions has 2nk + 4 paramters of which there are n — 1 
global gauge rotations, or 2nk — n + 5 gauge invariant paramters. 

Explicit solutions were found for the case of axial symmetry, an important limita- 
tion being the difficulty of solving the Nahm equation, or equivalently the quadratic 
ADHM constraint. We certainly expect more progress can be made on this in the 
near future. Nevertheless, we already found a rich structure that bodes well for being 
able to consider the constituent monopoles as independent objects. This is surprisingly 
subtle, as we have illustrated by the fact that an approximate superposition of charge 
1 calorons tends to give rise to a visible Dirac string. This is also related to the dif- 
ficulty of finding approximate multi-monopole solutions, which can be obtained from 
the caloron solutions by sending a subset of the constituent monopoles to infinity, as 
has been well studied in the charge 1 case ^ 

An important tool has been our study of the far-field limit, describing the abelian 
gauge field far removed from any of the constituent monopoles. This allows for a 
description of the long distance properties in terms of (self-dual) Dirac monopoles. 
Much could be extracted concerning its properties without the need to explicitly solve 
the Nahm equation. We conjecture in general to be able to identify the Dirac monopole 
location, but some work remains to be done here. 




Figure 6: The logarithm of the action density (cutoff for log(S') below -3) for an SU{2) 
charge 2 caloron with one type of constituent three times more massive than the other 
(/X2 = l/8, ai=eo" = 0,p, =2,e = 3.5) 

It may seem that all these results are somewhat academic since until recently none 
of these constituent monopoles were found in dynamical lattice configurations. First 
of all one would be tempted to search for them at high temperature, but it should be 
noted that above the deconfining phase transition the average Polyakov loop takes on 
trivial values, associated to the center of the gauge group, which is not the environ- 
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ment in which a caloron will reveal its constituents. This would give the well-known 



Harrington- Shepard solution constructed long ago |]14[. With our present understand- 



ing this solution can be seen as having n — 1 massless constituents which cannot be 
localized. Only when sending all these to infinity one is left with a monopole |jl5|. 
Furthermore, at high temperature classical configurations will be heavily suppressed 
due to their Boltzmann weight. Rather, the hope is that the constituent monopoles 
play an important role below the deconfining temperature, where the average of the 
Polyakov loop is non-trivial, and tends to favor equal mass constituents. This is why 
in this paper our examples were for that case, see Figs. |^J1,| and the discussion in 
Sect. p.2| . Nevertheless, the formalism developed here gives results for any choice of the 
holonomy, and a sample of unequal mass constituents is given in Fig. In general the 
constituent monopoles can be characterized by their magnetic (=electric) charge. For 
SU{n) there are n different types of abelian charges involved and all k constituents 
of a given type have the same mass. 

We will discuss briefiy the lattice evidence for the presence of constituent monopoles 
that has accumulated the last few years. A first numerical study using cooling was 
performed with twisted boundary conditions, which implies non-trivial holonomy ||T6 |. 
Good agreement was found with the infinite volume charge 1 analytic results, in par- 
ticular so for the fermion zero-modes |T^ which are more localized than the action 



density. A charge 2 solution was also found, shown in Fig. 8 of Ref. |]16l. Fermion 
zero-modes played an intricate role in an extensive numerical study of Nahm dualities 
on the torus . 



As suggested in Ref. [0] one may also enforce non-trivial holonomy on the lattice 
by putting at the spatial boundary of the box all links in the time direction to the 
same constant value Uq, such that U^^ = Voo (Nt the number of lattice sites in the 
time direction). This has been implemented in SU{2) lattice Monte Carlo studies as 
well, where Voo was set to the average value of the Polyakov loop, appropriate for the 
temperature at which the simulations were performed [^, Cooling was applied 
to find calorons, including those at higher charge. Apart from the configurations that 
in the continuum would be exactly self-dual, the lattice allows one to also consider 
configurations in which both self-dual and anti-selfdual lumps appear. This revealed 
constituent monopoles that seem not directly associated to calorons, called DD (as 
opposed to DD). Both objects in such a DD configuration have fractional topolog- 
ical charged, but opposite in sign. Perhaps these arise when two near constituent 
monopoles, one belonging to a caloron, the other to an anti-caloron, "annihilate" . Our 
analytic methods can not directly address this situation due to the lack of self-duality. 
The same holds for configurations that seem to only carry magnetic fields, which were 
already seen long ago P^ . 

One point of criticism that applies to both methods is that the choice of bound- 
ary conditions may force the "dissociation" of instantons into constituent monopoles. 
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particularly since volumes can not yet be chosen so large that many instantons are 
contained within a given configuration. A recent study [^] has done away with the 
fixed boundary conditions that enforce the non-trivial holonomy. Nevertheless, still 
one finds in many cases that calorons "dissociate" into constituent monopoles below 
the deconfinement transition temperature. A particularly useful tool has turned out 
to be the fermionic (near) zero-modes to detect the monopole constituents when they 
are too close together to reveal themselves from the action density ^3|. This re- 
lies on the observation that the zero-mode is localized on only one of the constituent 
monopoles, determined by the boundary conditions imposed on the fermions in the 
time direction [1^, |1^]. For SU{2) this is particularly simple, with periodic and anti- 
periodic boundary conditions of the fermions making the zero-mode switch from one 
constituent to the other, as is illustrated for a close pair of constituents in Fig. |^. In 
addition one may use the Polyakov loop for diagnostic purposes |2l|, for SU{2) 
taking the values I2 and — 12 near the respective constituent locations (at these points 
the gauge symmetry is restored, providing an alternative definition for the center of a 
constituent monopole). 



Figure 7: Fermion zero-mode densities as a function of t and z for a charge 1 caloron 
(/^2 = \ and p = \) with periodic (right) and anti-periodic (left) boundary conditions, 
compared to the action density (middle) (see also Fig. 5 in Ref. [^). Based on Ref. PB| ; 
produced with Ref. p^). 



Fermion eigenfunctions with eigenvalues near zero have also been used as an al- 
ternative to cooling, to filter out the high frequency modes and identify topological 
lumps. For a recent lattice study, including some discussion of calorons with non- 
trivial holonomy, see Ref. and references therein. Using the near zero-modes as a 
filter, constituent monopoles have even been identified recently for S'f/(3) well below 
the deconfining temperature , resembling Fig. |^ (see also Fig. 14 of Ref. 



For the exact axially symmetric multi-caloron solutions constructed in this paper 
the k associated fermion zero-modes (for charge k) will be derived in the near fu- 
ture. We anticipate that one can choose a basis where each is localized on one of 
the constituent monopoles, the type of which is determined by the choice of fermionic 
boundary conditions in the time direction. Analyzing these zero-modes is particularly 
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interesting in the light of some puzzles that were presented in a recent study |^ of the 
normalizable fermion zero-modes in the background of a collection of so-called bipoles, 
i.e. pairs of oppositely charged (but self-dual) Dirac monopoles, which are of interest in 
a wider context as well. This will be one of the many topics we have access to with our 
analytic tools. But ultimately our main aim is to develop a reliable method to describe 
the long distance features of non-abelian gauge theories in terms of monopole con- 
stituents to understand both confinement and chiral symmetry breaking. The results 
of this paper, and in particular the recent lattice results, provide some encouragement 
in this direction. 
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